A view of the narrow-line region in the infrared: active galactic 
nuclei with resolved fine-structure lines in the Spitzer archive 



K. M. Dasyra^'2, L. C. Ho^, H. Netzer^ F. Combes^, B. Trakhtenbrot^ E. Sturm^ L. 

Armus^, D. Elbaz^ 



O 
(N 

-)— » 
o 

O 

o 

(N 

o 
u 

6 

> 

a^ 
m 
m 

o 



% 



ABSTRACT 

We queried the Spitzer archive for high-resolution observations with the Infrared Spectro- 
graph of optically selected active galactic nuclei (AGN) for the purpose of identifying sources 
with resolved fine-structure lines that would enable studies of the narrow-line region (NLR) at 
mid-infrared wavelengths. By combining 298 Spitzer spectra with 6 Infrared Space Observatory 
spectra, we present kinematic information of the NLR for 81 z<0.3 AGN. We used the [Ne V], 
[O IV], [Ne III], and [S IV] lines, whose fluxes correlate well with each other, to probe gas pho- 
toionized by the AGN. We found that the widths of the lines are, on average, increasing with 
the ionization potential of the species that emit them. No correlation of the line width with 
the critical density of the corresponding transition was found. The velocity dispersion of the 
gas, cr, is systematically higher than that of the stars, tr*, in the AGN host galaxy, and it scales 
with the mass of the central black hole, Mbh ■ Further correlations between the line widths and 
luminosities L, and between L and Mbh, are suggestive of a three dimensional plane connect- 
ing log(MBH) to a linear combination of log(cr) and log(L). Such a plane can be understood 
within the context of gas motions that are driven by AGN feedback mechanisms, or virialized 
gas motions with a power-law dependence of the NLR radius on the AGN luminosity. The A/bh 
estimates obtained for 35 type 2 AGN from this plane are consistent with those obtained from 
the Mbh — cr* relation. 



Subject headings: galaxies: active — galaxies: kinematics and dynamics ■ 
Seyferts — infrared: galaxies — quasars: emission lines 



galaxies: nuclei — galaxies: 



1. Introduction 

Meaningful statistical tests of whether the bulk 
of black hole (BH) growth precedes, is paral- 
lel, or follows the peak of star-formation activ- 
ity need to rely on the comparison of several ob- 
servable parameters as a function of look-back 
time. In addition to the comparison of the star- 
form ation rate (SFR) with the BH accretion rat e 



(IR) bright galaxies. At low z and at low (Seyfert- 
like) lumi nosities, L 50-70% of the AGN are 
of tvpe 2 (IHo et al. |[l997bl : ISchmitt et al. |[200lt 
iHao et al. 1120051) . At z~I and at bolometric lumi- 
nosities that exceed 10^^ Lgun, at least I out of ev- 



ery 3 AGN is thought t o be of type 2 (jLacv et al. 



20071 : iGilh et aL 1120071 ) 



e.g., iMarconi et al. 1 12004 iMerloni et al. 1 120041) . 



the comparison of the stellar mass with the BH 
mass, Mbh, that is already accumulated at any 
given redshift z is also desirable. A significant 
impediment in performing the latter comparison 
comes from the fact that rapid BH growth often 
occurs in highly obscured environments, such as 
type 2 active galactic nuclei (AGN) and infrared 



In type 2 systems, the determination of A/bh 
is very challenging. Direct methods of measuring 
Mbh, e.g., by spatially resolv ed sub-parsec or par - 



sec scale kin ematics of stars ( Genzel et al. 19971) 



ionized gas ([Harms et al. Ill994i iMacchetto et al. . 
Il997l ). or water masers (JMivoshi et al. Ill995[ ). can 



only be applicable to very local objects. An al- 
ternative method that is success fully applied to 



type I AGN even at high z (.Kaspi et al. I I2007D 



uses the kinematics and the radial extent of gas 



clouds that are close enough to the BH to be 
gravitationally influenced by it. These clouds con- 
stitute the AGN broad line region (BLR). Emis- 
sion lines tracing the BLR, such as H/3, Mg II 
and C IV, have typical full width at half max- 
ima (FWHM) exceeding 2000 kms"!. The BLR 
radius is determined by reverberation mapping ex- 
periments, which use the natural variability of the 
AGN continuum and the time-delayed response of 
the gas (de)excitation to compute the light travel 
time acro s s the BLR (iBlandford fc McKee Ill982 : 
Petersoni Il993t iKaspi et al. I I2OOOI) . Assuming 
that the clouds are virialized, their distance from 
the BH and their velocity dispersion can be used to 
determine an enclosed dynamical mass and, thus, 
Mbh- However, the method is not applicable to 
type 2 AGN where the BLR is obscured. 

The lines that are seen in both type 1 and 
type 2 AGN are those that typically originate 
from gas clouds at a few p c to a few hundreds 
of pc away from the BH ( Netzer I l2004t iLaor 



20071 and references therein). These narrow-line- 
region (NLR) clouds have FWHM that are typ- 
icall y <1000 kms'^ Using the [QUI] 5007 A 
line, iNelson fc Whittle I ( 19961 ) demonstrated that 
the bulk of the gas traced by [O III] has a ve- 
locity dispersion that is comparable to that of 
the stars in the host-galaxy bulge. Subsequently, 
Mbh was foun d to scale with the width of the 
[O III] emission dNelson I I2OO0I: IShields et allbOOat 



Greene fc Ho 1 120051 ). A simila r result wa s found 
for the [N II] line at 6583 A byEo] llooi). 

Op tical NLR lines can suf fer severe obscura- 
tion ( Kauffmann et al. 1 120031) . This makes IR 
wavelengths more reliable and often unique for 
investigating relations between NLR kinematics 
and Mbh • To optimize the use of narrow lines as 
gravitational potential tracers for obscured AGN, 
we performed a similar analysis in mid-infrared 
(MIR) wavelengths using Spitzer and Infrared 
Space Observatory (ISO) high-resolution spectra. 
We demonstrated that the widths of the NLR lines 
[Ne V] at 14.32 i im. and [O IV] at 25 .89 /zm also 



scale with Mbh ( Dasvra et al. II2OO8I ). This result 



could provide a means of probing NLR kinematics 
and weighing BH masses in obscured galaxies at 
high z with the next generation IR telescopes. 

The goal of this paper is to perform an extended 
study of the NLR kinematics of optically selected 
type 1 and type 2 AGN in the MIR using all spec- 



tra that are available in the Spitzer archive. We 
aim to investigate for differences in the gas kine- 
matics as traced by various fine-structure lines, to 
test how the gas velocities compare with A/bh, 
and to further estimate the masses of local ob- 
scured AGN. For all computations, we use a 
AGDM cosmology with i/o=70 km s~^ Mpc~\ 
fi„i=0.3, and Oa=0.7. 

2. The Sample 

To study the NLR gas kinematics as probed by 
MIR emission lines, we queried the entire Spitzer 
archive for high-resolution observations of AGN 
obtained with the Infrared Spectrograph (IRS; 
Houck et al. 1120041 ). These spectra have a resolv- 
ing power of '^500 kms~^, sufficient for resolv- 
i ng features in the NLR of several local AGN 



(JDasvra et al. I I2OO8I ). We downloaded the re- 



served observations catalog (ROC) after the com- 
pletion of Spitzer' a cryogenic mission to ensure 
that all archival data that can be used for this 
study are included in our sample. We found 
1366 astronomical observation requests (AORs) 
containing high-resolution IRS spectra of extra- 
galactic sources performed in either single-target 
or cluster mode. 

We matched all IRS targets with optical spec- 
troscopic catalogs of AGN to identify type 1 
sources with existing A/bh estimates and [O III] 
5007 A detections that are required for the com- 
parison of NLR kinematics at various wave- 
lengths. For this purpose, we used i) the r e- 



verber a tion mapping catal o gs of iPeterson et al. 



(120041). iBentz et al 



s catalc 
T(l2006h 



and iDennev et al. I (|2010[ ). ii) the spectroscopic 



Bentz et al. 



(|2009f ). 



subsample of the Palomar-Green (PG) QSOs 



(IBoroson fc Green ll992HVestergaard fc Peterson 



20 Q6D, iii) the optica l spectr oscopic catalog o f 
Maraiani et al. I (|2003[ l. iv) the lHo fc Kim] (jiooi) 
catalog of type 1 AGN, v) the Sloan digital sky 
survey (SPSS) data re leas e 4 AGN catalogs o f 
Kauffmann et al. I (l2003l). oflGreene fc Ho I (l2005l). 
and of iNetzer fc Traktenbrot I (|2007t ). and vii) 
the SDSS data release 7 (DR7) . To identify 
AGN in the DR7, we used the [Baldwin et al. I 
(|l98ll ) diagnostic diagram of [O III]/H^ ver- 
sus l^N II] /Hg with the AGN boundaries as up- 



dated by iKewlev et al. I (|2006l) . We also exam- 



ined sources that have a QSO target type, or 



a broad-line or AGN spectral type subclass as- 
signed by the SDSS pipeline. These samples 
also include a few type 2 AGN, which we further 
co mplemented by matchin g all IRS targe t s wit h 



thelMulchaev et al. 



Bassani et al. 
Reyes et al 



19941) .iTurner et al. 



(2008) 
2009|) catalogs 
optical spectra of t ype 2 



(■1999^ 'Zakamska et aL 
Bennert et al 



Liu et al. 



(20091) 



1997 ). 



20031) . 



Spinelli et all (|2006l ) 



A 
AGN 



few 



and 
additional 
found in 



In total, we found and reduced datasets from 
370 Spitzer AOKs, corresponding to high-resolution 
spectra of 298 sources. Resolved NLR lines were 
determined from the inspection of the [S IV] 
10.51 Aim, [Nelll] 15.56 /um, [O IV] 25.89 /im, 
or [Ne V] 14.32 /im profiles as described in Sec- 
tion [321 To these 298 AGN, we added 6 sources 
with fine structure lines that were similarly re- 
solved by IS O short wavelengt h spectrograph 



(SWS) data (jSturm et al. I 120021: iDasvra et al. 



I2OO8I ). The redshift distribution of ah 304 sources 
in the final sample is shown in Figure [TJ The 42 
type 1 AGN and the 39 type 2 AGN that have re- 
solved MIR narrow lines are presented in Tables [1] 
andlU respectively. 

3. Data Reduction 

3.1. Spectral Extraction 

We downloaded the Spitzer S15.0 pipeline basic 
calibration data (BCD) files. We used both staring 
and mapping observations, taken on either single- 
source or cluster mode. 

The BCD files of the long-high (LH) wavelength 
data of each source were processed by the IDL 
routine DARK SETTLE, which is posted at the 
Spitzer Science Center (SSC) web site, to correct 
for gradations of the dark current along the detec- 
tor that lead to order tilting and mismatch. Using 
the short-high (SH) BCD files and the dark-settle 
corrected LH BCD files, we computed the aver- 
age frame for each set of on-source observations, 
namely, for each each target, nod position, and 
module. To identify cosmic ray hits, we also com- 
puted the median frame for the same set of ob- 
servations and compared it to the average frame 
on a pixel by pixel basis. A pixel was flagged if 
the difference between the median and the average 
frame value differed by a factor of more than three 
times the root-mean square (rms) noise in the me- 
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Fig. 1. — Redshift distribution of AGN with MIR [S IV], 
[Ne III], [O IV], or [Ne V] detections (grey solid hne), and 
redshift distribution of AGN with resolved [S IV], [Ne III], 
[O IV], or [Nc V] lines (black dotted line). 

dian frame. For this pixel, the average value was 
replaced by the median value. For sources with 
available observations on sk}i], we also computed 
the median sky frame, which we subtracted from 
the on-source frame. 

The next step was the removal of bad and rogue 
(i.e., slowly varying time response) pixels. We 
used the individual sky frames to create a generic 
bad pixel mask for all sources observed in a sin- 
gle Spitzer campaign, which we merged with the 
bad pixel map available at the SSC web site for 
the same campaign. We then merged this generic 
mask with the mask of each individual (on-source 
and on-sky) exposure to create the mask for each 
science frame, nod position, and module. We pro- 
ceeded to further masking of outliers, i.e., pixels 
that were located a couple of columns away from 
the edges of each spectral order, and that had a 
value exceeding the rms noise of the science frame. 
Their values were replaced with the median value 
of the frame, computed using only pixels in the 
useful detector range. A final visual inspection 
and manual cleaning of the science frames was per- 
formed using the IRSCLEAN routine. All rogue 
pixels and outliers were fiagged in the mask file. 



^ Several programs whose scientific goals required only the 
use of line fluxes (but not the use of equivalent widths or 
continuum flux measurements) did not acquire sky obser- 
vations given that the flux of a line does not depend on its 
underlying continuum level. 



The uncertainty of the average on-source frame 
was calculated as the square root of the sum of the 
squares of all individual uncertainty files, divided 
by the number of exposures. When sky observa- 
tions were available, we computed the sky frame 
uncertainty in a similar way and combined it with 
the on-source frame uncertainty to produce the 
uncertainty of the final science frame. 

We used the science frame, together with its un- 
certainty and mask files as input to the SSC soft- 
ware SPICE, which produces a one-dimensional 
spectrum from a two-dimensional spectral image. 
To extract the spectra, we used the regular extrac- 
tion mode, which equally weighs pixels when col- 
lapsing them along each row. We assumed a point 
source extraction, since the LH and SH apertures 
(22'.'3xll'.'l and ll'.'3x4'.'7, respectively) are likely 
to include the bulk of the NLR emission for most 
sources in our sample. The end product of SPICE 
is the wavelength and flux calibrated spectrum for 
each individual order. 

We merged the spectra of the various orders 
to a single spectrum for each nod position, clip- 
ping noisy edges (between 2-25 pixels, depending 
on the order). This task was performed for both 
the SH and the LH datasets. The SH and LH spec- 
tra were then merged to produce the full-range 
spectrum per nod position. The final spectrum of 
each object was produced by averaging the one- 
dimensional, full-range spectra of the two nod po- 
sitions. Special care was taken not to merge sky- 
subtracted and non-sky-subtracted datasets for 
sources with observations from different programs. 
At the wavelengths of key atomic/molecular lines, 
only nod positions without bad pixels were used, 
when possible. For all other pixels with a bad pixel 
flag in one nod position, we only kept the value of 
the second nod position when the two values dif- 
fered by more than the local rms noise value. 

3.2. Line Fitting 



As in iDasvra et al. I (|2008[ l'^l. we fitted all lines 
and their underlying continua using Gaussian 
and second-order polynomial functions, respec- 
tively. Gaussian fitting was preferred over direct 
FWHM computations which, given the resolu- 




Fig. 2. — Stacked profiles (dotted lines) constructed us- 
ing 29 2<0.02 individual sources with unresolved line de- 
tections at all wavelengths. Solid lines correspond to the 
Gaussian best-fit solutions of the stacked profiles. The dif- 
ferences in the line peak positions are indicative of wave- 
length calibration uncertainties. 



tion of IRS, can only be reliably performed for 
few high- velocity systems (e.g., ISpoon fc Holt 



2009). Gaussian fitting also suits best the studies 
of cloud kinematics on virial, regular, or symmet- 
ric motions, instead of clouds that are entrained 



by asymmetric outflows (Gre ene fc Ho Il2005f) 



'AH previously presented data were reanalyzed to account 
for changes in the Spitzer pipeline, and for enhancements 
in our bad-pixel identification routines. 



To claim a line detection, we required that its 
signal-to-noise ratio (S/N) was greater than 3. To 
measure line widths, however, we only used lines 
with S/N>5 to avoid studying profiles of barely 
detected lines. We considered resolved all lines 
with CTm — £„! > cn + ^i, where am is the measured 
velocity dispersion, ai is the instrumental resolu- 
tion at a given wavelength divided by 2.35, and e^ 
is the error of (Ti . The average resolution value in 
the 12.0—18.0 /im range, which comprises all neon 
lines for the low z galaxies, is 507±66 kms"^. 
The 25.0—34.2 /xm range, which comprises most 
of the O IV emission, has an average resolution of 
503±63 kms~^. The error of a^, Cm, encapsu- 
lates both measurement and instrumental uncer- 
tainties. It was computed as (e^j -I- ef)'^'^, where 
est is the standard deviation of the different veloc- 
ity dispersion values that were obtained for each 
line when using different polynomial functions to 
describe its underlying continuum. Intrinsic veloc- 
ity dispersions, a, were computed as (tr^ — erf )°'^, 
converted to restframe, and presented in Tables [1] 
and [2] for key fine-structure lines in all sources. 

To assess possible systematic errors in the wave- 
length calibration, we normalized and stacked the 
line profiles of 29 sources with detected but un- 
resolved emission lines of [Ne II] , [S IV] , [Ne III] , 
[O IV], and [Ne V] (Figure [2]). The stacking of un- 
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Fig. 3. — Relation of the [Ne V] luminosity, eonimonly used as an AGN indicator, with the luminosities of other fine-structure 
lines in MIR and optical wavelengths. The optical luminosities are not corrected for extinction. 



resolved lines in different-z sources permitted for a 
finer sampling of the instrumental resolution curve 
than that determined by the detector pixel size. 
We only used sources at z<0.02 to study system- 
atic effects within a few pixels from the restframe 
wavelength of each line. This z cutoff translates to 
25 pixels at the center of the SH array range, 14.7 
/im. We found that line peak offsets due to such 
systematics are typically limited to 50 kms^^. 

4. Results 

4.1. Identifying NLR Tracers in the MIR 

Of the 304 AGN in our sample, 300 had spec- 
tral coverage at 10.51 ^m, and 143 showed [S IV] 
emission. Similarly, 226 sources had spectral cov- 
erage at 25.89 /xm, and 135 had [O IV] detections. 
Thus, the typical line detection rate was of order 
-50%. 

We plotted the luminosities of the most fre- 
quently observed MIR lines against that of [Ne V] 
to identify fine-structure lines that can be reliably 
used as tracers of the gas that is photoionized by 
the AGN (Figure [Sj. Given that photons of at 
least 97.12 eV are required to ionize Ne IV to Ne V, 



its ionization source must be far-ultraviolet or soft 
X-ray radiation from the AGN. We find that the 
[Ne V] luminosity is on a tight correlation of 0.2 
dex scatter wi th the [O IV] and [S IV] lum inosi- 
ties (Figure [51 iPereira-Santaella et al. Il2010l ). and 
we further confirm that the [Ne V] and [Ne III 



luminosities trace each other well ( Goriian et al. 



20071 ). Therefore, the Ne V, O IV, Ne III, and S IV 



ions are most likely to be primarily excited by the 
same mechanism. However, the scatter between 
the luminosities of [Ne V] and [Ne II] or [Si II] is 
approximately twice as high, of 0.4— 0.5 dex. Ne II 
and Si II have low ionization potentials (21.56 and 
8.15 eV, respectively), thus, a non- negligible part 
of their line emission could be tracing star- forming 
regions in the AGN host galaxy. 

The excellent, almost one-to-one correlation 
over several orders of magnitude in luminosity that 
is seen between the [Ne V] 14.32 /xm and [S IV] 
10.51 /xm lines does not include any extinction 
correction. This result indicates that the bulk 
of the [S IV] emission is not strongly affected by 
the silicate absorption feature centered at 9.7 ^m, 
which is seen in several of our targets with mod- 
erate or high optical depths, i.e. r>0.5. Thus, 



the geometric distribution of the sihcates is hkely 
to be m ore compact than the size of the NLR 
(see also'Soifer et al. 1 120021: iTristram et al. II2007I: 
ISchweitzer et al. 2008! ) . While obscuration effects 
do not significantly affect the relative fluxes of 
MIR lines, the [O III] 5007 A emission can suffer 
from strong extinction primarily in type 2 AGN 
(Figure IS]). 

4.2. Line Profiles: Velocity Dispersion of 
Different Ionized Gas Components 

Our query for sources with resolved [S IV], 
[Ne III], [O IV], or [Ne V] lines resulted in 81 AGN 
with kinematic information of the NLR gas in 
the MIR (Figures H and O. The different fine- 
structure lines that were used as primary tracers 
of clouds photoionized by the AGN often led to 
different cr measurements (Figures S] and [S]) . Dif- 
ferences in the ionization potential of the various 
ionic species can play a major role in determin- 
ing the observed line profiles, alongside with dif- 
ferences in the critical density for collisional de- 
excitation of the various transitions, and with light 
extinction by dust particles. We find the [Ne V] 
line to be often broader than lines from ions of 
lower ionization potential (e.g., Mrk590; Figures |4] 
and m Tables [3] and |4]) . A significant increase in 
the average velocity dispersion with increasing ion- 
ization potential is shown in Figures [5] and [71 as 
found for 16 sources with resolved profiles in all 
[Nell], [OIII], [SIV], [Nelll], [O IV], and [Ne V] 
lines. This result provides evidence for a stratifi- 
cation of the NLR clouds, i.e., for ions of high ion- 
ization being located near er to the BH than ions 



of low ionization potential (IFilippenko fc Halpern 



1984 lOliva et al. 111994 iHo 1120091 ). 



counter example to this statistical finding could be 
MCG-03-34-064. Because all of its neon lines have 
higher velocity dispersions than [O IV] and [S IV] , 
its NLR density could be of order 10* hydrogen 
atoms per cubic centimeter. 

Dust either mixed or outside the NLR gas 
clouds can lead to different extinction of lines 



Radially dependent line widths will also be pro- 
duced if the NLR has a gas dens ity gradient in- 

creas i ng with proximity to th e BH (|Filippenko fc Halpern 



1984 [Ferguson et al. 1 119971 ). which is not con 



nected to a mean ionization level gradient. The 
line widths will be affected as the line fluxes from 
transitions of low critical density for collisional 
de-excitation, ric, will be preferentially suppressed 
at dense environments. We observe no trend of 
(7 with Tic (Figure [6]) , suggesting that the typical 
NLR gas densities are below ric for the MIR tran- 
sitions that we examined. The lack of dependence 
of the average a on ric demonstrates the power of 
MIR lines in probing the NLR kinematics. One 



at di fferent A ([Groves et al. 112004 [Winter et al. 
20101 ). Because MIR lines are less affected by ob- 



scuration than optical lines, they could be probing 
clouds that are nearer to the BH, where obscu- 
ration is often higher. In this case, their widths 
could be broader than those of optical lines, in 
particular in type 2 AGN. Mrk273, NGC2623, 
and IRAS15001+1433 do have larger [O IV] than 
[O III] widths and low [O III]/ [O IV] flux ratios 
(see Figure [31), but their [O III] and [S IV] or 
[Ne III] widths are in good agreement with each 
other. The comparable ionization potentials of 
O III, S IV, and Ne III, i.e., 35.12, 34.79, and 
40.96 eV respectively, ascribe again any profile dif- 
ferences to ionization effects. We conclude that 
obscuration is not a common driver of differences 
in the profiles of optical and IR lines. The nu- 
clear obscuration could be high enough to hide 
blue wings, associated with outflows moving away 
from the observer, even in IR wavelengths. 

A simple comparison of the NLR gas velocity 
dispersion to that of the stars in the host galaxy 
shows that a is not identical to tr* (Figure ^. 
The gas velocity dispersion systematically exceeds 
that of the stars, with a large scatter between the 
two quantities. The average velocity dispersion 
excess is in the range 50-100 kms~^ for all lines. 
The excess could be intrinsically higher, given that 
the optical observations had a slit width that was 
often an order of magnitude narrower than the 
IRS slit width. The origin of this excess could 
be related to either gas clouds that axe acceler- 
ated by AGN- feedback mechanisms ( Ho [ [20091 ) , 



or to virialized gas clouds with a more compact 
spatial distribution than that of the stars, when 
probed by emission lines from ions of ionization 
potential > 35eV. High spatial resolution MIR 
spectroscopy of 8 local AGN indicated that the 
[S l y] emission is unresolved at a scale of 100 



pc ([Honiget al. [[2008[ ). while NIR integral-field- 



unit data of Circinus showed that most of the 
[Si VI] and [Ca VIII] ernission is unresolved at 4 pc 



(JMiiller Sanchez et al. [[20061 ). This velocity dis- 
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Fig. 4. — Profiles of fine-structure lines in type 1 AGN with at least one resolved NLR line in the MIR. For each source, 
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Figure (4]- continued. 
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Fig. 5. 



Profiles of fine-structure lines in type 2 AGN with at least one NLR line resolved in the MIR. 
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Figure [5]- continued. 
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Fig. 7. — Stacked MIR line profiles of the 16 sources 
used in Figure |6] 
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Fig. 6. — Average NLR gas velocity dispersion (corrected 
for instrumental broadening) as a function of ionization po- 
tential (upper panel) and critical density in units of elec- 
trons Tie per cubic centimeter (lower panel) for 16 sources 
with detections in all of the lines used for the construction 
of this figure. 



persion excess is seen even for the [O III] 5007 A 
line, which is resolved for all systems in our sam- 



ple. iGreene &: Ho 



(|200a ) had also reported such 
an excess in a large flux-limited sample of AGN 
observed with SDSS. 

Another trend known from optical NLR studies 
is that the widths of the lines increase with their 



own luminosities () Phillips et al. I Il983t [Whittle 



1985lll992b[) We reproduce this result for the MIR 
fine-structure lines that we examined in Figure [5] 
The velocity dispersion correlates roughly as L*^'^^ 
based on, e.g., the [Ne V] measurements. A higher 
dependence of tr on L is plausible when taking into 
account the upper limits of a. An exact determi- 



nation of this slope will require the use of a flux- 
limited sample of AGN in the local universe. 



4.3. Line Profiles: 

tries 



Shifts and Asymme- 



Changes in the profile moments of different 
lines can also be seen when the emission orig- 
inates from gas out of dynamical equilibrium. 
Several of the 81 sources in our sample show 
signatures of outflowing (or inflowing) gas mo- 
tions, such as an offset from systemic velocity. 
This offset often depends on the ionization po- 
tential of each ionic species. With increasing ion- 
ization potential, the lines trace clouds that are 
nearer to the BH, and therefore more suscepti- 
ble to outflows. We identified 6 sources in which 
the offset between the peak of the lowest ioniza- 
tion line available and the highest ionization line 
available exceeds Nyquist sampling of the resolu- 
tion element, >250kms^^. Namely, these are the 
type 1 AGN 3C273, IRAS13342-H3932, and the 
type 2 AGN IRAS05189-2524, IRAS15 001-H433, 
IRAS23060+0505, and Mrk6 09 (see also lSpoon et al. 
20091 : ISpoon k HoiTI l2009l ). Line shifts below 
Nyquist sampling of the instrumental resolution, 
yet above the wavelength calibration uncertainty 
(50 kms~^), could also be theoretically observed. 
They could be considered real if they systemat- 
ically increased as a function of the ionization 
potential of the ionic species. Only one sourcqj 



'^Other sources, such as NGC2623, NGC7469, Arpl48, and 
PG1211-)-143, could also have outflows of ionized gas, which 
need to be proven with spectroscopy of higher resolution 
than that of IRS. The blue wings of [Ne V] and [O IV] in 
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satisfies this criterion, Mrkl457. Its [Ne V] and 
[Ne II] lines are offset by 100 kms^^. 

Asymmetric wings are found in '^l/5th of 
the sources with resolved profiles. We consider 
such wings to be reliable only when they are 
detected in two or more lines. Similar [O III] 
5007 A and [Ne III] or [Ne V] wings are observed, 
for example, in PG1351+64G, PG1411+442, and 
PG1440+356. Likewise, an agreement of the 
[O III] and [O IV] 25.89 /.tm line profiles is found 
for 3C273, IRAS13451+1232, IRAS15001+1433, 
MrklG14, NGC4235, and NGC7674. Specifically, 
IRAS1345H-1232 is a merging syst em with two 
nucle i separated by ^5 k pc (e.g., lAxon et al. I 
I2OOOI : iDasvraet aUbOOeal) . The velocity of the 
secondary [O IV] peak , at ~-1000kms~^, is com- 
parable to the velocities of the outflowing [O I] 
and [O III] gas components that are seen in op- 
tical spectroscopy and that are associated with 
the nucleus respo nsible for the radio-jet emission 



The slope was set to 4.24, which is identical to 
that of the stellar relation (Ferrarese fc Merritt 



()Holt et al. Il2003h . 



5. Discussion 

5.1. A Three Dimensional Plane connect- 
ing Mbh? and the L and a of the NLR? 

By performing a similar profile analysis of the 
[Ne V] and [O IV] lines, restricted only to AGN 
with direct A/bh measurements from r everbera- 
tion experiments (jPeterson et al. 1 120041) . we pre- 



viously demonstrated that the NLR velocity dis- 
persi on correlates with the mass of its central 
BH (JDasvra et al. 1 120081 ). We now further pop- 



ulate this relation using data from the full Spitzer 
archive, complementing them with ISO SWS data 
to cover the parameter space a< 200 km s^^. The 
Mbu values for this expanded sample include 
single-epoch optical spectroscopic estimates for 
type 1 AGN, as well as direct Mbh measurements 
for a few type 2 AGN (see Section [21 Tables [3] 
andg]). 

The fit to these data (Figure \W\i was per- 
formed with the IDL routine MPFITFUN. Given 
the sparsity of data at the low a end, we opted 
for a linear fit of fixed slope in logarithic space. 



NGC6240 are also redshifted with respect to [Ne II] and 
[Ne III] , but their overall line shifts are hard to determine 
as these lines are b lended with [CI II] and [Fe II] (see also 
lArmus et"ani2006t) . 



I2OOO I: 'Gebh ardt et al. "2000*) as recently revisited 
bv iGiiltekin et al. I ( 2009) . While the rms scatter 
was computed using all available datapoints, the 
best-fit solution was computed using only Spitzer 

A/bh>10^M(7) or 



^ 



observations of sources with 
ISO observations, as in iDasvra et al. I ()2008[ ). 
The reason why we excluded Spitzer datasets 
for AfBH<lO^M0 is that they would introduce 
a bias toward high intercept values: the reso- 
lution of IRS is insufficient to resolve lines on 
the left-hand side of the A^bh — cr relation be- 
low this threshold. This excludes the type 2 AGN 
NGC1068, and the narrow-line Seyfert Is (NLSls) 
]VIrk493, and NGC4051 that have high a values 
for their BH masses. The existence of such out- 
liers nonetheless suggests that a could fail as a 
proxy of A/bHi as it is also known for the BLR 



gas fe.g.. lVestergaard fc Peterson 1120061 ). On the 



other hand, such sources are not outliers in the 
relation connecting A^bh to the luminosity of the 
NLRhnes (Figure [IJ). 

The multiple relations between A/bh and a, a 
and L, and A/bh and L could be suggestive of 
a plane connecting these three parameters (Fig- 
ure [12]), which we fit using the equation 

log{ Mbh) = alog{a) + lilog[L) + 7, (1) 

where a,/3, and 7 are constants. We find that the 
best-fit-solution coefficients correspond to a=0.9 
and /3=0.5, when averaged over all IVIIR lines. 
Given the L and (T~^ proportionality shown in Fig- 
ure [HI this result roughly reproduces the A/bh ~ 
a^ relation. In Figure [131 we present the plane 
equation for each line, when fixing the log(i7) and 
log(-L) slopes to their average values for simplic- 
itjo- We find that the use of a plane equation 
minimizes the scatter of the NLR-based A/bh es- 
timates from their actual values. The improve- 
ment primarily originates from the correction of 
the A/bh — o" relation outliers at its low- A/bh and 
high- cr end. 



^^ We have adapted this approach throughout all fitting pro- 
cedures in this work (see also Figures fTol and fTTjl to facilitate 
comparisons of the rms scatter among the various relations, 
and to avoid biases related to the small galaxy sample that 
is used for their creation and that varies from line to line. 
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5.2. The Origin of the Scaling Relations 
between the NLR Gas Properties and 

Mbh 

The physical interpretation of the suggested 
plane linking the NLR line luminosity, veloc- 
ity dispersion and the BH mass, varies depend- 
ing upon the assumed kinematics and distribu- 
tion of the gas clouds. A plane equation can 
be meaningful for clouds on accelerated mo- 
tions powere d directly by the AGN , via r adiation 



pressure (Murayama fc Taniguchi I 119981 ) acting 



mos tly upon dust partic les due to their high opac- 



ity (iDopita et al. 1 120021) . magnetic fields related 



to jets (|Whittlelll992cD . and AGN-related winds. 



The winds can lead to either asymmetric outflows, 
identified in 7% of the systems in this sample, or 
to symmetric outflow s that can be related t o the 
AGN accretion disk ( Crenshaw et al. I l2003t IHo 
20091 ). In this scenario, the measured gas veloc 



ity would be related to the fraction of the energy 
generated by the AGN that is deposited into the 
NLR gas. The luminosity would be a probe of 
the AGN accretion rate, since the luminosities of 
MIR NLR lines are known to be correlated with 
the opti cal. X-ray, and bolome tric AGN luminosi- 
ties (S chweitzer et al. i 2006; Dasvra et al. I l2008t 



Melendez et al. IbOOSHRigbv et al II2009I : R. Mor 
et al. in preparation). A plane connecting a, 
L, and Mbh would then suggest that the kinetic 
energy of the interstellar medium (ISM) , as mea- 
sured from the MIR lines, is directly related to the 
amount of material that is accreted onto the BH 
for a given Mrh value. An analo gous idea was 



introduced bv lMerloni et al. I ( 20031 ). who found a 



correlation between the mass of an accreting BH 
with its X-ray and radio luminosities. 

Instead of responding to AGN feedback mech- 
anisms, the gas clouds could be in virial motions 
that are dictated by the enclosed mass at the NLR 
radius R. For a geometric factor / converting this 
total mass to the BH mass, the virial equation can 
be written as 

log{ A/bh) = 2logia) + log{R) + log{f) - Zog(G), (2) 

where G is the gravitational constant. In this 
model, the NLR radius depends on the AGN lu- 
minosity alongside with the stellar mass and dis- 
tribution in the AGN host galaxy. Equation ^ 
will take the form of Equation ([IJ if i? scales with 
the AGN luminosity as a power law of I/, as found 



for the [O III 



.5007 A line (JBennert et al. 112002 : 

Schmitt et al. 1 120031 ) . The hardness of the ioniz- 



ing radiation together with the ionization fraction, 
the covering factor, and the density distribution 
of the clouds can also determine how far the AGN 
radiation reaches, illuminates, and photoionizes 
gas clouds. For example, the observed increase 
of a with L can be driven either by a tendency 
for m ore massive BHs to reside in larger galaxies 
(JHo I 2009) , or by the ionization state of the NLR. 
Lines from ions of high ionization potential are 
thought to be tracing matter-bound clouds that 
are partially ionized (iMuravama fc Taniguchi I 






19981 : IWilson et al. 1119971 ) . In the case of partially 
ionized clouds, an increase of the AGN luminosity 
can lead to an overall expansion of the NLR, while 
most of the ionization can still occur for clouds at 
small radii. The exact relation between a, L, and 
Mbh will depend on the gas density distribution, 
which is encapsulated in the scaling factor /. The 
coefficients of the best-fit plane solution, a=0.9 
and /3=0.5, suggest that / is dropping with cr, 
corresponding to a lower scaling factor for denser 
gas within a fixed radius. If however, the scaling 
factor / was constant, a would be equal to 2, and 
the best-fit plane solution would correspond to 
/3=0.4. Such a solution could also be plausible. It 
would increase the rms scatter by a small amount, 
i.e., by only 0.01-0.02 dex for the MIR lines and 
by 0.07 dex for [O III]. 

If the cloud kinematics cannot be approximated 
by either feedback-driven motions or by virial mo- 
tions. Equation [1] might not provide a good means 
of estimating A/bh- For cloud kinematics that are 
described by a combination of virial and acceler- 
ated motions, e.g., in a virialized NLR where radi- 
ation pressure is significant, Mbh is proportional 
to a^RG^^ + a'L, where a' is a constant that de- 
pends on the c olumn density and leve l of ionization 
of the clouds (iMarconi et al. I l2008t but see also 



Netzer fc Marziani 1 120101) . Further investigation 
of the optimal description of Mbh from narrow 
line properties will be tested on large, flux-limited 
samples of AGN with multiwavelength datasets. 
This will be the focus of a forthcoming paper us- 
ing SDSS data. 
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Fig. 8. — Comparison of the stellar veloeity dispersion with the NLR gas velocity dispersion, as measured from MIR lines 
and from the optical [O III] 5007 A line. 
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5.3. Black Hole Mass Estimates in Ob- 
scured AGN 

Given that we see no significant difference in 
the NLR properties of ty pe 1 and type 2 AGN at 
a giv en L or a^, (see also iPereira-Santaella et al. 



20101 ) , we applied the relations presented earlier in 
this work to type 2 AGN for which such an anal- 
ysis is presently possible from Spitzer data. Only 
four of these sources have other, direct A/bh rnea- 
suremcnts (Table [2]) • The comparison between the 
Mbh estimates from the stellar Mbh — cr* rela- 
tion, the A/bh — cr relation for the NLR gas, and 
the best-fit linear combination of log(cr) and log(L) 
are presented in Table [5] For either computation 
based on the NLR gas, the result is averaged over 
all available MIR lines. We find that the me- 
dian ratio of A/bh as estimated from the plane 
equation over A/bh as estimated from the stel- 
lar velocity dispersion is 1.8. On the other hand, 
the median ratio of A/bh as estimated from the 
NLR gas velocity dispersion over A/bh as esti- 
mated from the stellar velocity dispersion is much 
higher, 6.5. The difference seen when folding a lu- 
minosity dependence on the computation of A/bh 
is largest for the IR-bright galaxies IRAS05189- 
2524, IRAS15001-I-1433, and Mrk609, which have 
gas motions that are predominantly out of dynam- 
ical equilibrium. 

6. Conclusions 

We queried the full Spitzer archive for high- 
resolution IRS spectra of type 1 AGN with BLR- 
based A/bh estimates from optical spectroscopy, 
and type 2 AGN. We analyzed the spectra of 298 
objects, which we combined with ISO spectra of 
6 more AGN, aiming to study the fundamental 
properties of the NLR in the IR, and to calibrate 
its fine-structure-line widths and luminosities to 
the mass of the central black hole. We found the 
following. 

• All of the [Ne V], [O IV], [Ne III], and [S IV] 
line luminosities are well correlated with 
each other (with an rms scatter of ~0.2 dex). 
Given the high ionization potential of Ne IV 
to Ne V, their emission is attributed to (fur- 
ther) photoionization of ions by the AGN. 
Thus, all of these lines can be used to study 
the NLR kinematics. 



• The bulk of the silicates responsible for the 
9.7 /xm absorption feature is likely to be lo- 
cated in a region smaller than the NLR. Oth- 
erwise, the luminosity of [S IV] at 10.51 /im 
would not correlate equally well with that of 
all other IR narrow lines without an extinc- 
tion correction in obscured AGN. 

• Of the 304 sources in our sample, 81 had at 
least one [NeV], [O IV], [Ne III], or [S IV] 
narrow line that was resolved. We find that, 
on average, the line widths increase with in- 
creasing ionization potential of the species 
that emit them in the range 22— 97eV. Dif- 
ferent fine-structure lines probe different 
(locations within the) clouds, with high- 
ionization-potential ions being preferentially 
found nearer to the BH. 

• No trend was observed between the line 
width and critical density for these transi- 
tions, indicating that the average density of 
the NLR gas is typically below 10* hydrogen 
atoms per cubic centimeter. 

• The velocity dispersions of all MIR lines that 
we examined, as well as that of the [O III] 
5007 A line, are systematically above the ve- 
locity dispersion of the stars in the AGN host 
galaxies. Moreover, the line widths increase 
with their own luminosities, as known from 
optical wavelengths. 

• The mass of the central BH correlates with 
the velocity dispersion and with the lumi- 
nosity of the NLR lines. We speculate that 
these results could be suggestive of a three- 
dimensional plane, connecting log(A/BH) to 
a linear combination of log(cr) and log{L). 
Such a plane equation could be meaning- 
ful for a virial distribution in which the 
NLR radius has a power-law dependence on 
the AGN luminosity, or for AGN-feedback 
driven motions. 

• Given that we find no significant differences 
in the NLR kinematic properties of type 1 
and type 2 AGN, we used the NLR gas prop- 
erties as seen in the MIR to estimate the BH 
masses residing in 35 local obscured AGN. 
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Table 1 
Black hole masses & optical NLR and host galaxy properties of type-1 AGNs. 
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87 ± 5 


NGC4235 


0.0080 






39.53 




170 


143 ± 34 


NGC5548 


0.0172 


8.20 ± 0.00 




41.44 




210 


192 ± 15 


NGC7469 


0.0165 


1.22 ± 0.14 




41.54 




153 


142 ± 3 


PG0026 + 129 


0.1455 


39.30 ± 9.60 




42.46 




191 




PG0804 + 761 


0.1010 


69.30 ± 8.30 




42.06 




393 




PG1119 + 120 


0.0502 


2.95 ± 0.59 




41.49 




249 


162 ± 28 


PG1229+204 


0.0630 


7.32 ± 3.52 




41.68 




197 


162 ± 32 


PG1351+640 


0.0882 


36.30 ± 7.26 




42.19 




390 




PG1411+442 


0.0896 


44.30 ± 14.60 




41.82 




393 




PG1426+015 


0.0866 


129.80 ± 38.50 




41.89 




287 


217 ± 15 


PG1440 + 356 


0.0790 


2.94 ± 0.59 




41.48 




393 




PG1613 + 658 


0.1295 


27.90 ± 12.90 




42.08 




354 




PG2130+099 


0.0630 


45.70 ± 5.50 




41.93 




296 


172 ± 46 


PG2349-014 


0.1742 


196.00 ± 19.60 




43.11 




221 


223 ± 36 


3C120* 


0.0330 


5.55 ± 2.70 




41.62 




122 


162 ± 24 


Mrk509* 


0.0344 


14.30 ± 1.20 




42.33 




221 




NGC37S3* 


0.0097 


2.98 ± 0.54 




41.43 




98 


95 ± 10 


NGC4151* 


0.0033 


4.57 ± 0.52 




41.44 




181 


95 ± 7 


Note.— (*) sourcos with MIR f 


ne-structure lines resolved o 
measurements are taken fro 


nly by 

Othe 
ted fr 
spectr 


ISO. 
rberation mappi 
rwise, thoy arc 
om single-epoch 


ng oxpcrimc 

taken from 

spectrosco 






Col. (3): When possible, A^rm 


nts IPe 


erson et al. 1120041; 








poch spectroscopy 








pv IMa 


zianietal.l 120031; 


INetzeriTraktenbrot 1120071: IHo 


fe Kim 1120091: SDSS archival 


al as inlVester« 


ard & Pete 


rsoni l200ell. 


2MASXJ10614428 + 3539304 has v 


™ak, bro 
|lM2aD, 
resented 


ad /fa detection b 


rt no 


H^ detection. 
Hot Klml nroSD. SDSS 


DR7, as 




Col. i4): Data from 1 Whittle 1 


IMarziani et al.| f 


well as on 


. The 


measurements from the spectra p 




J|199 


^, performed as 


inlNetzer fe Trakte 


fT^T^nnnnn 


luminosities are not corrected for 


extincti 
r fromE 
spersion 

] PTOTTD 


an. 






sured from 


SDSS SF 




Col. (5): ^[o III] ™1™" ™ t-ke 


hittiei iiiSis;iai).rTTT 

and the 7th data 


fc Kim 1 120091), or mea 

aken fromlWelson et al 

release of SDSS. Wh 


ectra. 


Dasyra 
its were 


Col. (6): Host galaxy velocity di 


]|2(w;^),;^ 


nken et 


iT] \-iini/\\ 


et al. (2006a; 2006b; 2007), [He 


n more tha 


n one n 




available, their average value v/ 


as used. 
3 combin 


IRAS13451 + 1232 
3d value for both ni. 


cloi. 


erger lUasyra ei 


al. llSijijfial 


. The 


velocity dispersion 


presented here corresponds to th 
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Table 2 
Black hole masses & optical NLR and host galaxy properties of type-2 AGNs. 



Galaxy 


z 


'l^fBH 


log(L[Qjjjj/(crgs s-1)) o-[Q jjj,(5007A) 


^. 


(-) 


(-) 


(10^ Mq) 


<-) 


(km=-l) 


(km.-l) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


2MASXJ08035923+2345201 


0.0297 




40.99 


121 


130 ± 5 


2MASXJ0S244333+2959238 


0.0254 








41.13 


157 


107 ± 5 


2MASXJ101S192S+3722419 


0.0497 








41.33 


155 


100 ± 6 


2MASX J 12384342 + 0927362 


0.0829 








41.97 


233 


229 ± 10 


2MASXJ16164729+3716209 


0.1518 








42.51 


150 


229 ± 9 


CGCG218-007 


0.0273 








40.85 


131 


149 ± 5 


ESO103-G035 


0.0133 














IRAS05189-2524 


0.0430 








41.82 


230 


137 ± 16 


IRAS15001 + 1433 


0.1623 








41.52 


193 


296 ± 19 


IRAS18325-5926 


0.0202 














IRAS23060 + 0505 


0.1730 








42.24 


243 




MCG-03-34-064 


0.0165 














Mrkl066 


0.0120 








40.88 


172 




Mrkl457 


0.0486 








41.38 


182 


143 ± 6 


Mrk273 


0.0378 








40.79 


255 


244 ± 18 


Mrk3 


0.0135 








42.14 


362 




Mrk463E 


0.0510 








42.64 


140 


163 ± 8 


Mrk609 


0.0345 








41.32 


215 


145 ± 4 


Mrk622 


0.0232 








40.48 


291 


135 ± 5 


NGC1068 


0.0038 


0.86 


+ 


0.03 


41.80 


451 


151 ± 7 


NGC1275 


0.0174 








41.68 


601 


246 


NGC2622 


0.0286 








41.43 


255 




NGC2623 


0.0185 








39.01 


140 


152 ± 9 


NGC2639'' 


0.0111 








39.70 


226 


188 


NGC3079'' 


0.0037 








37.72 


527 


150 


NGC425S 


0.0018 


3.78 


+ 


0.01 


40.87 


174 


167 


NGC4507 


0.0118 








41.53 


91 




NGC5256 


0.0279 








41.91 


176 


187 ± 8 


NGC5506 


0.0062 








40.58 


106 




NGC572S 


0.0094 








41.12 


136 




NGC5929 


0.0083 








40.02 


179 


123 ± 3 


NGC6240 


0.0245 








39.02 


426 


229 ± 43 


NGC7172 


0.0087 














NGC7674 


0.0290 








42.01 


189 




SBS1133 + 572 


0.0516 








41.47 


168 


207 ± 7 


UGC02608 


0.0233 














UGC5101 


0.0394 








39.91 


193 


188 ± 6 



0.0018 
0.0014 



4.90 
0.17 



± 1.40 
± 0.03 



150 ± 7 
SO 



Note. — (*) : 
(**) sources v 
Col. (3): Mhh 



ourcc's with shallow optical spcctr; 
ith MIR fine-structui 



nd |N¥ 



ertainty. 
Col. (4): Th. 



et al. I (^TO7|). 



a lLodato' 



ir spectra obtained u 

only by ISO. 

i Berlin I 120031 1. R?T^ 



nhill et aTH 120031 ^. fllei 
available, we used the 



In 



i iUJflTaH .rns 



Col. (5): <T[o jij] vain 
SDSS spectra. 
Col. (6); Host galaxy 
IJasvra et aT~| |iiOOab| l 



sn 



puted from IWiiittle | ^^^) , [UT 
d Ir 



s. 



et al. I U554I) , Veilleu 
a. They are no 
et al. ItrgSai . lWliittle'l |ltit>2a|l . [V511eux et al. 1 (155311 . 



. et ; 



I] I2005II . 



|ii(Wab| l. rFacon-tiarrOBo-gn 
:V I liiMl . maiteUnetal.i l 



piled from ll'alirTI | 



.[Te 



■ith the 
(1995; 1999), 
-ed from 
et al. 



,EI] 

d from the 7th data release of SDSS. 



ran nT!T!7j ). n'gESoni et ai. i rot 

'..,'l'. .' .'l I I.' I |M 11??^ 



24 



Table 3 
Fluxes and Velocity Dispersions of Ionized Gas Narrow Lines in type 1 AGNs. 













































Galaxy 


i'?" "I , 


If '''' 2 
10 ^° W m~^) 


f[Nc III] 
(10^^^ W m^^) 


fjO IV] 


IL"'^ V] 


"[Ne 


",1 


•^[S IV] 


"[N 


e III] 


''[O IV] 


"jNe 


XI 




(-) 


(10~-^° W m~^) 


(10^^® W m^^) 


(10 ^^ \V m ^) 


(km s 


-1) 


(km 


s 


-1) 


(km 


s" 


-1) 


(km 


= -1) 


(km 


s" 


-1) 




2MASSJ09184900+2117170 


< 3.61 


< 10.86 


15.65 ± 2.39 




< 25.07 












206 


± 


35 














2MASSiJ1659397+183436 


< 10.67 


< 10.16 


26.79 ± 2.95 




< 10.15 












277 


± 


52 














2MASXJ 105 14428+3539304 


< 6.21 


9.97 ± 1.98 


19.45 ± 2.11 




< 5.20 






325 


± 


48 




















3C249.1 


4.37 ± 0.40 


8.85 ± 0.43 






5.05 ± 0.59 


370 ± 65 


286 


± 


33 












329 


± 


36 




3C273 


11.42 ± 1.67 


36.04 ± 3.33 


41.58 ± 2.61 


85.45 ± 6.42 


21.20 ± 4.21 


230 ± 39 


470 


± 


87 


402 


± 


47 


450 


± 67 


495 


± 


57 




Arkl20 


28.97 ± 2.53 


< 13.56 


35.23 ± 2.71 


36.84 ± 3.70 


10.04 ± 1.61 


202 ± 32 








213 


± 


38 


190 


± 39 










CGCG121-075 


24.88 ± 1.44 


24.39 ± 2.00 


47.28 ± 1.96 


85.25 ± 1.86 


16.36 ± 1.25 






209 


± 


50 


238 


± 


42 


220 


± 27 


223 


± 


43 




BSO140-G043 


96.32 ± 2.92 


87.41 ± 2.95 


139.5 ± 2.2 


247.6 ± 4.9 


79.10 ± 2.19 






















203 


± 


29 




FAIRALL9 


21.50 ± 1.85 


28.59 ± 3.86 


44.80 ± 3.44 


63.73 ± 3.68 


20.79 ± 3.77 






296 


± 


45 


197 


± 


36 


258 


± 48 


231 


± 


42 




IC4329A 


241.5 ± 6.2 


245.7 ± 9.5 


535.0 ± 6.8 


1009 ± 14 


286.1 ± 8.4 


227 ± 32 


235 


± 


29 


256 


± 


30 


275 


± 31 


282 


± 


34 




IRAS13342 + 3932 


53.68 ± 1.11 


19.74 ± 1.27 


< 28.42 


98.54 ± 1.72 


33.98 ± 0.98 






310 


± 


48 








324 


± 29 


403 


± 


38 




IRAS13451+1232 


45.66 ± 1.82 


< 5.59 


< 70.31 


27.87 ± 5.24 


< 7.86 


311 ± 37 














479 


± 54 










IRAS1S216 + 6418 


27.52 ± 2.06 


68.26 ± 1.23 


104.5 ± 3.5 


240.1 ± 12.3 


50.79 ± 1.69 


214 ± 34 


304 


± 


31 


237 


± 


31 


343 


± 36 


390 


± 


34 




MCG- 2-58-22 


74.71 ± 1.86 


< 33.31 


87.08 ± 3.79 


126.3 ± 3.6 


28.26 ± 2.63 












183 


± 


33 


187 


± 32 


204 


± 


43 




Mrkl014 


58.26 ± 1.42 


33.20 ± 1.67 


93.20 ± 2.96 


119.9 ± 5.8 


46.01 ± 2.39 


253 ± 37 


376 


± 


36 


379 


± 


41 


379 


± 39 


380 


± 


48 




Mrk279 


81.42 ± 2.88 


26.48 ± 2.04 


80.35 ± 3.97 


102.8 ± 3.6 


34.80 ± 2.19 


215 ± 34 


237 


± 


36 


222 


± 


35 






248 


± 


42 




Mrk335 


10.54 ± 0.90 


< 25.80 


22.56 ± 1.75 


130.0 ± 3.7 


12.88 ± 1.86 


















155 


± 17 


247 


± 


64 




Mrk493 


62.69 ± 1.31 


< 23.47 


27.09 ± 1.53 


28.42 ± 4.76 


9.59 ± 1.27 
















280 


± 


39 


267 


± 46 


280 


± 


38 


to 


Mrk590 


32.81 ± 2.12 


< 16.58 


21.14 ± 2.54 


31.04 ± 3.06 


9.57 ± 1.64 
















197 


± 


33 


200 


± 33 


206 


± 


40 


Mrk704 


< 10.05 


< 61.44 


< 50.82 


< 151.67 


51.08 ± 3.04 


























249 


± 


37 




Mrk705 


52.30 ± 1.90 


< 23.26 


49.59 ± 2.11 


57.82 ± 2.19 


30.18 ± 1.66 
















192 


± 


31 


210 


± 35 


280 


± 


39 




NGC3227 


708.7 ± 25.2 


220.9 ± 5.5 


725.8 ± 6.2 


668.0 ± 17.6 


231.2 ± 7.5 










241 


± 


32 


228 


± 


30 


216 


± 29 


209 


± 


34 




NGC3516 


74.01 ± 2.29 


128.9 ± 4.4 


164.4 ± 2.9 


458.4 ± 5.5 


67.57 ± 2.18 
















198 


± 


30 


229 


± 28 










NGC4051 


172.6 ± 4.9 


< 47.59 


162.0 ± 2.9 


347.6 ± 7.1 


107.7 ± 3.9 
















172 


± 


29 


295 


± 27 










NGC4235 


34.85 ± 1.62 


< 6.03 


33.29 ± 1.77 


33.60 ± 3.04 


< 5.28 
















214 


± 


34 


286 


± 43 










NGC5548 


83.08 ± 3.32 


42.90 ± 2.55 


81.76 ± 2.77 


124.7 ± 8.6 


31.50 ± 2.16 
















177 


± 


30 






206 


± 


38 




NGC7469 


1915 ± 27 


90.00 ± 7.88 


357.9 ± 7.5 


340.0 ± 38.0 


154.4 ± 10.0 










158 


± 


12 


196 


± 


29 


196 


± 45 


240 


± 


44 




PG0026-I-129 


2.56 ± 0.33 


4.94 ± 0.43 


< 7.72 


20.16 ± 4.01 


< 5.07 






















305 


± 40 










PG0804-I-761 


< 4.44 


15.53 ± 2.28 


20.89 ± 1.70 


20.74 ± 2.87 


< 11.13 










292 


± 


111 


323 


± 


47 


342 


± 66 










PG1119-t-120 


< 4.00 


20.63 ± 2.55 


28.41 ± 2.01 


60.00 ± 2.37 


16.16 ± 2.33 


























193 


± 


47 




PG1229-I-204 


5.84 ± 0.59 


< 20.90 


12.05 ± 1.37 


26.74 ± 3.55 


11.48 ± 1.13 


























218 


± 


SO 




PG1351-t-640 


21.92 ± 0.99 


< 10.26 


29.39 ± 1.73 


< 12.38 


< 3.83 
















200 


± 


43 














PG1411-I-442 


4.54 ± 0.79 


7.79 ± 1.00 


9.24 ± 0.62 


13.99 ± 2.31 


6.49 ± 0.58 






















212 


± 45 


316 


± 


63 




PG1426-t-015 


13.46 ± 1.08 


13.67 ± 2.28 


25.55 ± 1.09 


29.51 ± 4.38 


8.47 ± 1.09 


226 ± 38 








258 


± 


36 


292 


± 57 


289 


± 


57 




PG1440-I-356 


42.72 ± 0.96 


15.55 ± 1.70 


38.57 ± 1.26 


50.22 ± 2.37 


14.26 ± 0.86 


193 ± 34 


211 


± 


39 


253 


± 


35 


253 


± 38 


304 


± 


107 




PG1613-t-658 


37.66 ± 1.23 


13.25 ± 1.01 


30.90 ± 1.02 


66.49 ± 3.15 


< 7.92 


226 ± 30 


269 


± 


33 


289 


± 


33 


225 


± 29 










PG2130-I-099 


16.14 ± 1.22 


32.54 ± 1.53 


55.73 ± 2.88 


100.8 ± 3.7 


45.16 ± 3.01 












187 


± 


30 






287 


± 


45 




PG2349-014 


15.95 ± 0.71 


< 6.36 


20.51 ± 1.01 


34.73 ± 3.82 


< 6.83 


273 ± 48 








259 


±_ 


34 


308 


± 35 










3C120 


78.86 ± 3.72 


225.7 ± 3.3 


271.6 ± 3.7 


1129 ± 7 


163.8 ± 2.7 


















114 


± 8 


111 


T 


16 




Mrk509 


118.0 ± 4.1 


73.19 ± 3.92 


153.0 ± 3.4 


180.0 ± 5.1 


54.04 ± 2.63 


173 ± 35 














165 


± 30 










NGC3783 


197.0 ± 3.3 


128.8 ± 5.1 


255.9 ± 5.0 


380.0 ± 8.5 


152.0 ± 5.2 












244 


± 


12 


141 


± 24 


133 


± 


21 




NGC4151 


1180 ± 12 


1130 ± 22 


350.0 ± 13.9 


2030 ± 33 


560.0 ± 11.0 






207 


± 


16 


155 


± 


7 


169 


± 12 


139 


± 


25 




Note— No data tor the lino flux 


cs indieate no speetra 


coverage at the 




lengths of those lines. 


No data tor the veloc 


ty dispc 




of dote 


ct 


jd lir 


res ind 


eato that the 


r signal 


to no 




ratio 




was between 3 and 5, or that the 


ir profiles were unreso 


Ived. 







































Table 4 
Fluxes and Velocity Dispersions of Ionized Gas Narrow Lines in type 2 AGNs. 



C-) 



^[Ne II] 



*'[S IV] 
-18 w ^ 



*'[Nc III] 



f[0 IV] 

(10-18 T^ m-2j 



*'[Nc V] 



^[Ne II] 



^[S IVJ 



^[No III] 
km s-^) 



•^[O IV] 
(km s-1) 



^[Nc V] 



to 



2MASXJ08035923+2345201 

2MASXJ08244333+2959238 

2MASXJ 10 18 1928 + 3722419 

2MASX J 12384342 + 0927362 

2MASXJ 16 164729 + 3716209 

CGCG218-007 

ESO103-G035 

IRAS05189-2524 

IRAS15001 + 1433 

IRAS1S325-5926 

IRAS23060 + 0505 

MCG-03-34-064 

Mrkioee 

Mrkl457 

Mrk273 

Mrk3 

Mrk463E 

Mrk609 

Mrk622 
NGC1068 
NGC1275 
NGC2622 
NGC2623 
NGC2639 
NGC3079 
NGC4258 
NGC4507 
NGC5256 
NGC5506 
NGC5728 
NGC5929 
NGC6240 
NGC7172 
NGC7674 
SBS1133 + 572 
UGC02608 

UGC5101 

Ccntaurus A 

CircinuB 



< 6.39 


< 5.50 


18.78 ± 2.08 


28.45 ± 1.59 


< 4.75 














203 


± 


42 












33.17 ± 2.04 


< 20.88 


49.47 ± 3.44 


83.64 ± 2.13 


36.14 ± 2.78 


207 


± 


39 




















264 


± 30 


8.67 ± 1.53 


< 6.35 


14.69 ± 1.87 


28.81 ± 1.60 


16.59 ± 2.80 


























200 


± 32 


< 8.33 


20.09 ± 2.25 


27.85 ± 2.15 


58.83 ± 2.38 


12.94 ± 2.35 








232 


± 


42 


260 


± 


42 


261 


± 


36 


265 


± 73 


< 5.34 


16.53 ± 1.49 


22.55 ± 2.52 


90.14 ± 2.28 


17.88 ± 2.23 


























211 


± 37 


102.6 ± 3.1 


35.72 ± 2.64 


86.07 ± 2.69 


166.7 ± 2.3 


46.59 ± 1.98 


























215 


± 33 


294.0 ± 7.6 


115.1 ± 10.1 


414.0 ± 9.4 


311.7 ± 13.9 


172.9 ± 10.5 


223 


± 


37 


232 


± 


34 


276 


± 


31 


240 


± 


36 


346 


± 52 


191.8 ± 8.9 


65.30 ± 6.26 


186.1 ± 7.1 


261.2 ± 37.2 


152.6 ± 12.0 


253 


± 


42 


335 


± 


55 


392 


± 


41 


444 


± 


63 


404 ± 43 


66.48 ± 1.36 


4.74 ± 0.82 


27.97 ± 0.69 


21.78 ± 3.21 


15.89 ± 1.97 


247 


± 


34 


388 


± 


36 


348 


± 


38 


558 


± 


96 


422 


± 81 


367.9 ± 7.8 


114.9 ± 8.2 


436.0 ± 5.3 


390.7 ± 12.2 


262.6 ± 9.6 


181 


± 


31 


265 


± 


39 


367 


± 


32 


300 


± 


34 


404 


± 64 


32.26 ± 2.40 


19.65 ± 2.34 


25.61 ± 2.10 


36.31 ± 6.45 


20.12 ± 2.69 


278 


± 


59 


366 


± 


127 


285 


± 


62 


343 


± 


67 


419 


± 62 


514.9 ± 8.2 


476.1 ± 8.3 


1150 ± 13 


1053 ± 22 


616.6 ± 8.2 


344 


± 


34 


298 


± 


33 


361 


± 


28 


273 


± 


30 


435 


± 36 


1094 ± 21 


102.0 ± 7.4 


469.1 ± 7.6 


418.1 ± 22.1 


92.38 ± 8.31 














216 


± 


32 


284 


± 


41 


246 


± 48 


91.99 ± 2.39 


< 14.09 


38.13 ± 4.09 


28.10 ± 2.54 


22.09 ± 3.04 




















228 


± 


43 


310 


± 60 


444.9 ± 7.9 


101.7 ± 2.3 


338.1 ± 2.5 


560.1 ± 17.8 


112.3 ± 3.7 


252 


± 


30 


357 


± 


34 


332 


± 


30 


357 


± 


38 


406 


± 32 


979.8 ± 10.2 


592.7 ± 6.2 


1749 ± 11 


1964 ± 24 


632.5 ± 7.5 


276 


± 


29 


322 


± 


29 


314 


± 


29 


263 


± 


28 


327 


± 31 


108.2 ± 3.5 


275.0 ± 5.4 


404.6 ± 7.3 


641.8 ± 11.3 


193.9 ± 4.7 


308 


± 


35 


284 


± 


30 


236 


± 


29 


296 


± 


29 


287 


± 31 


213.6 ± 6.2 


< 22.18 


56.71 ± 2.53 


82.43 ± 26.69 


38.69 ± 4.28 














202 


± 


36 








488 


± 60 


97.29 ± 13.80 


< 8.61 


49.29 ± 2.82 


< 37.42 


< 38.06 














333 


± 


34 












4988 ± 183 


6199 ± 336 


13781 ± 246 


20406 ± 473 


8974 ± 221 


359 


± 


52 


414 


± 


44 


383 


± 


31 


364 


± 


48 


363 


± 39 


461.6 ± 8.0 


< 16.26 


223.7 ± 5.6 


< 84.37 


< 11.66 


299 


± 


31 








260 


± 


31 












62.02 ± 2.63 


< 18.46 


80.41 ± 2.93 


89.29 ± 3.58 


24.73 ± 2.08 














213 


± 


33 








289 


± 39 


550.1 ± 9.7 


11.20 ± 1.36 


147.2 ± 1.7 


117.4 ± 18.2 


37.55 ± 2.97 














213 


± 


31 


357 


± 


52 


222 


± 39 


86.23 ± 2.66 


< 6.90 


48.88 ± 1.95 


< 23.82 


< 7.51 


346 


± 


43 








336 


± 


49 












1310 ± 60 


< 14.61 


238.3 ± 3.3 


< 94.19 


< 38.79 


198 


± 


33 








285 


± 


29 












123.5 ± 7.3 


< 15.61 


70.61 ± 5.61 


79.40 ± 8.03 


< 13.67 


193 


± 


42 








191 


± 


36 


204 


± 


53 






307.8 ± 6.9 


93.08 ± 6.84 


287.8 ± 6.3 


354.1 ± 11.6 


125.7 ± 6.4 














220 


± 


35 


318 


± 


42 


345 


± 81 


160.6 ± 3.1 


27.06 ± 1.56 


94.40 ± 1.45 


569.8 ± 10.3 


20.54 ± 1.37 


227 


± 


29 


239 


± 


37 


227 


± 


30 


191 


± 


28 


202 


± 40 


850.6 ± 14.3 


735.4 ± 15.6 


1537 ± 11 


2262 ± 40 


568.0 ± 10.6 








239 


± 


33 


207 


± 


28 








208 


± 32 


366.2 ± 7.7 


317.2 ± 8.2 


536.5 ± 6.7 


1155 ± 12 


217.1 ± 3.4 








272 


± 


41 


263 


± 


29 


306 


± 


30 


314 


± 32 


104.8 ± 2.2 


12.41 ± 1.98 


95.07 ± 1.40 


43.67 ± 6.18 


< 6.09 














210 


± 


28 












1872 ± 29 


38.7 ± 2.9 


599.9 ± 4.2 


374.0 ± 35.9 


90.24 ± 9.02 


362 


± 


37 


279 


± 


52 


341 


± 


29 












334.6 ± 8.1 


55.74 ± 4.43 


168.1 ± 4.0 


384.4 ± 5.3 


90.23 ± 3.49 


256 


± 


30 








204 


± 


30 












215.7 ± 5.2 


152.6 ± 5.9 


346.0 ± 6.5 


442.9 ± 12.1 


184.4 ± 7.0 








267 


± 


34 


253 


± 


32 


242 


± 


32 


267 


± 39 


108.3 ± 2.9 


42.01 ± 3.62 


105.3 ± 3.8 


168.7 ± 2.8 


69.89 ± 2.87 








248 


± 


40 


















566.5 ± 9.4 


292.8 ± 9.5 


706.3 ± 6.1 


1378 ± 17 


312.9 ± 5.1 








213 


± 


31 


















379.7 ± 9.0 


< 7.34 


146.2 ± 2.1 


57.18 ± 9.66 


34.19 ± 3.73 


254 


± 


40 








350 


± 


34 


285 


± 


41 


385 


± 40 



2210 ± 45 
4536 ± 145 



140.0 ± 9.6 
1270 ± 58 



1400 
4000 



± 12 
± 90 



980.0 ± 42.2 
6793 ± 420 



200.0 ± 8.5 
2180 ± 101 



134 ± 21 
86 ± 4 



127 ± 9 
113 ± 3 



132 ± 21 
90 ± 8 



Tablc[3] 



Table 5 
Estimates of A/bh in type-2 AGN. 



(10*^ Mq) 



(10' Mq) 



Ai"BH (plane equation) 
(10^ Mq) 



2MASXJ08035923 + 2345201 


2.1 


2MASXJ0S244333 + 2969238 


0.9 


2MASXJ101S1928 + 3722419 


0.7 


2MASXJ123S4342 + 0927362 


23 


2MASXJ 16 164729+3716209 


23 


CGCG21S-007 


3.8 


ESO103-G035 




IRAS05189-2524 


2.7 


IRAS15001 + 1433 


70 


IRAS18325-5926 




IRAS23060+0505 




MCG-03-34-064 




Mrkioee 




Mrkl457 


3.2 


Mrk273 


31 


Mrk3 




Mrk463B 


5.5 


Mrk609 


3.4 


Mrk622 


2.5 


NGC1275 


32 


NGC2622 




NGC2623 


4.1 


NGC2639 


10 


NGC3079 


3.9 


NGC4507 




NGC5256 


9.9 


NGC5506 




NGC572S 




NGC5929 


1.7 


NGC6240 


23 


NGC7172 




NGC7674 




SBS1133 + 572 


15 


UGC02608 




UGC5101 


10 



6.6 
7.2 
29 
132 



9.2 
8.8 



4.4 
7.5 
7.3 



7.3 
8.4 



7.9 
10 



9.0 
11 
8.2 
7.5 
4.7 
3.0 
8.4 
8.1 
6.0 
8.7 
3.0 



Note.— The i 



third ( 



alu 



in Figure^] The fourth 
gas- baaed estimate 
i\'/gp^ uncertainties 
higher, typically w 



Mgpj estin 



present the 
af a factor 



Late using the A'. 
the plane cquati' 



of the stellar 
BH - ^ relatit 



>f 0.5 dex. 
aagnitude. 



Mgpj — CT relation. The 
n for the NLR, as found 
Figure rT2l For the NLR 
ines. Formal, statistical 
: the uncertainty can be 
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